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Abstract. We develop a diagrammatic approach for solving few-body problems in heteronuclear fermionic 
mixtures near a narrow interspecies Feshbach resonance. We calculate s-, p-, and (i-wave phaseshifts for 
the scattering of an atom by a weakly-bound dimer. The fermionic statistics of atoms and the composite 
nature of the dimer lead to a strong angular momentum dependence of the atom-dimer interaction, which 
manifests itself in a peculiar interference of the scattered s- and ^waves. This effect strengthens with 
the mass ratio and is remarkably pronounced in ^°K-(^°K-^Li) atom-dimer collisions. We calculate the 
scattering length for two dimers formed near a narrow interspecies resonance. Finally, we discuss the 
collisional relaxation of the dimers to deeply bound states and evaluate the corresponding rate constant 
as a function of the detuning and collision energy. 
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1 Introduction 

Numerous advances in the field of ultracold Fermi gases 
over the past decade have enabled the exploration of novel 
strongly interacting regimes in fermionic systems (see [1] 
and [2] for review). The BCS-BEC crossover, extensively 
studied in either Potassium (^^K) or Lithium (^Li) ho- 
monuclear systems, is now being actively pursued in the 
new generation of experiments on mixtures of these two 
isotopes [3-10]. The mass ratio is thus a new parameter 
introduced into the crossover phase diagram [11]. It is then 
natural to ask whether changing this parameter can lead 
to qualitatively new crossover physics and, if so, how large 
should the mass ratio be in order to see non-trivial effects? 

For a very large mass ratio (of the order of several 
hundreds) a crystalline phase can emerge on the BEG side 
of the crossover [12]. The effect is due to a strong long- 
range repulsion between the heavy fermions originating 
from the exchange of their light partners. Another mass 
ratio dependent change in the behavior of the system oc- 
curs in the problem of two identical heavy fermions of 
mass interacting resonantly with a light atom of mass 
171^. The exchange of the light atom results in an attrac- 
tive potential between the heavy fermions proportional to 
1/m^i?^, where R is the distance between them. This ex- 
change attraction competes with the repulsive centrifugal 
barrier oc for the identical fermions. For mass ra- 

tios m^/mi larger than the critical value 13.6 the exchange 
attraction dominates over the centrifugal barrier and the 



heavy particles experience a fall to the centre in the 1 / R^- 
potential [13]. On resonance (interspecies scattering length 
a = oc) this leads to the Efimov effect - the existence of an 
infinite number of bound heavy- heavy- light trimer states 
[14-16]. For smaller mass ratios the centrifugal barrier is 
dominant. On the one hand this effective three-body re- 
pulsion excludes the Efimov effect. On the other it sup- 
presses recombination processes requiring three atoms to 
approach each other to very short distances, which is evi- 
dently advantageous for the collisional stability of the gas. 
The lower the mass ratio, the more stable this three-body 
system is [17]. 

It turns out that even for m^/m^ < 13.6 the tt>l- 
system exhibits non-perturbative effects on the positive 
(BEG) side of the resonance, where there is a weakly 
bound heteronuclear molecular state. One of us found 
that the three-body recombination to this state vanishes 
for m^/m^ ^ 8.6 [18]. Later, Kartavtsev and Malykh ar- 
gued that this phenomenon is related to the existence of 
a weakly bound not Efimovian trimer state for / > 
8.2 [19]. The trimer has unit angular momentum and for 
smaller mass ratios turns into a p-wave atom-dimer scat- 
tering resonance. We have recently shown that in the case 
of a K-Li mixture (m^/m^ = 6.64) the K-(K-Li) atom- 
dimer scattering should be dominated by this p-wave res- 
onance in a wide range of collision energies [20]. Moreover, 
by introducing an external quasi-2D confinement, the p- 
wave atom-dimer interaction can be tuned from attractive 
to repulsive, allowing for a trimer formation. 
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In this paper we develop a uniform-space diagram- 
matic approach for studying few-body processes in a hete- 
ronuclear fermionic mixture near an interspecies Feshbach 
resonance of finite width. We calculate relevant atom- 
dimer scattering phaseshifts and partial cross-sections in 
the homonuclear case and in the K-Li case. Passing from 
m^jm^ = 1 to m^jm^ = 6.64 we observe an increase in 
the atom-dimer interaction, repulsive in even angular mo- 
mentum channels and attractive in odd ones. The most 
dramatic increase is found in the channel with unit angu- 
lar momentum - the p-wave scattering volume changes by 
more than an order of magnitude. Our exact calculations 
are complemented by a qualitative explanation of the ob- 
served effect based on the Born-Oppenheimer approach, 
which we generalize to the case of a narrow interspecies 
resonance. We predict a very strong interference between 
s- and p-waves in atom-dimer scattering. Depending on 
the collision energy, the scattering is dominant in back- 
ward or forward directions, which can be observed exper- 
imentally by colliding an atomic cloud with a cloud of 
molecules. We use our diagrammatic approach to calcu- 
late the dimer-dimer scattering length a^d as a function 
of the atomic scattering length a and the width of the 
interspecies resonance. Finally, we discuss the main mech- 
anisms of the collisional relaxation of dimers into deep 
molecular states, and calculate the corresponding atom- 
dimer and dimer-dimer relaxation rate constants as func- 
tions of a, the width of the resonance, and the collision 
energy. 

The paper is organized as follows. In Sec. 2 we dis- 
cuss the two-body problem in the narrow resonance case 
and introduce our field- theoretical approach. The main 
part of the paper is structured according to the previous 
paragraph - Sees. 3 and 4 are devoted to the three- and 
four-body problems respectively. In Sec. 5 we discuss the 
inelastic collisional relaxation in atom-dimer and dimer- 
dimer collisions, and in Sec. 6 we conclude. 



2 Two-body problem near a narrow Feshbach 
resonance 

We assume that all interatomic interactions in the t-i 
fermionic mixture are characterized by van der Waals po- 
tentials. We also assume that the intraspecies interactions 
are not resonant, and therefore can be safely neglected 
in the ultracold regime. Let us denote the van der Waals 
range of the interspecies interaction by R^. and write down 
the partial wave expansion of the on-shell scattering am- 
plitude [13] 



/(k, kO = ^(2^ + l)P^(cos Q)h{^)- 



(1) 



^=0 



Here k and k' are initial and final relative momenta such 
that |k| = |k^| = /c, and d = Zk,k' is the scattering angle. 
We set /i = 1. The partial wave amplitudes /^(A:) can be 
written in terms of the phase shifts (^^(/c) as 



Expanding the denominator of Eq. (2) in powers of kR^ 
gives the effective range expansion. In particular, in the 
s-wave channel = 0) we have 



A:cot^o(^) ~ ^ + 2^^^^ 



(3) 



and the corresponding expansion in the j9-wave channel 
= 1) reads 

A:^cot(^i(A:) -v'^ -h ^^k^k^ + . . . , (4) 

where v is the p-wave scattering volume, and k^ is a pa- 
rameter analogous to the effective range. 

The scale of the scattering length a, the effective range 
ro, and other expansion parameters in the higher order 
terms in Eq. (3) are set by the length R^ or its power 
of suitable dimension, and the same holds for higher par- 
tial waves. For A: ^ the partial scattering amplitudes 
are proportional to {kRe)'^^. Thus, in the limit kRe <C 1 
(ultracold regime) the s-wave scattering amplitude, which 
equals /o(0) = —a, is the most important interaction pa- 
rameter in the mixture. 

Near a scattering resonance the scattering length can 
be modified and, in particular, can take anomalously large 
values (i.e. \a\ ^ R^). A magnetic Feshbach resonance oc- 
curs when the collision energy of the two atoms is close 
to the energy of a quasidiscrete molecular state in another 
hyperfine domain, which is called closed channel. The tun- 
ing of the scattering amplitude is achieved by shifting the 
open and closed channels with respect to each other in an 
external magnetic field (hyperfine states corresponding to 
the open and closed channel have different magnetic mo- 
ments). The width of the resonance is determined by the 
strength of the coupling between these two channels. The 
narrower the resonance, the stronger the collision energy 
dependence of the scattering amplitude, and, therefore, 
the larger the effective range ro- We call a resonance nar- 
row, if |ro| ^ Re [21]. In fact, near such a resonance ro 
is necessarily negative and it is convenient to use another 
length parameter [22] 



-ro/2 



1 



2/iabg/irel^^ ' 



(5) 



where /i = m^m^/ {m^ + m^) is the reduced mass, is 
the background scattering length, /i^ei is the difference in 
the magnetic moments of the closed and open channels, 
and AB is the magnetic width of the Feshbach resonance. 
All ^Li-^^K interspecies resonances discussed so far are 
characterized by i?* > lOOnm [4,7], which is much larger 
than the van der Waals range Re ~ 2.2nm. 

One can imagine an interatomic potential for which the 
higher order terms in Eq. (3) are also anomalously large. 
For example, we can introduce one or several additional 
closed channels with quasistationary states very close to 
the open-channel threshold resulting in a rather exotic 
scattering amplitude [24]. However, in this paper we as- 
sume a more practical and simple case in which the terms 
denoted by ... in Eq. (3) vanish in the limit kRe 0. 
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We will also assume that scattering with ^ > is not reso- 
nant and can be neglected in this limit. Then, substituting 
Eq. (3) into Eq. (2) we get the well-known formula for the 
resonant scattering at a quasidiscrete level [13] (written 
as a function of momentum rather than energy) 



(6) 



The ratio i?* /a measures the detuning from the resonance 
and we distinguish the regime of small detuning, i?*/a <C 
1, and the regime of intermediate detuning, i^*/a ^ 1. 
The properties of a few-body system in these two limits 
are qualitatively different [22]. 

In order to describe the mixture near a narrow res- 
onance we use the two-channel Hamiltonian [25] 



2m, 



-jf (&p«E+k,taE_k,; + &P«p-k,i«f +k,t) ' (^) 



k,cr=t4 



where a| ^ and are creation and annihilation operators 

of the two fermionic species while (b) creates (annihi- 
lates) a closed- channel molecule of mass M = + m^. 
The atom-molecule interconversion amplitude g is taken 
constant up to the momentum cut-off A oc and ujq 

is the bare detuning of the molecule. The quantities a and 
i?* are related to the parameters of the model (7) by [23] 
(see also Appendix A) 



1^9 



1 



(8) 



The bare propagators of atoms and closed- channel mo- 
lecules read 



Gt,^(p,Po) 
^o(p,Po) 



po -p2/2m^4 -hiO' 
1 

Po-p^/2M -ujo^iO' 



(9) 



where -\-iO slightly shifts the poles of G and Dq into the 
lower half of the complex po-plane. A physical dimer con- 
sists of a closed-channel molecule dressed by open-channel 
atoms. The corresponding propagator is given by (see Ap- 
pendix A) 



D{p,Po) 



27r//i 



.(10) 



-iO 



The pole of D(0,po) determines the dimer binding energy 

eo = -( VTT4R*7« - 1) VSm^*^- (11) 

Equation (11) interpolates between the two limits: for 
small detuning we have eg — — 1/2/ia^ and in the regime 
of intermediate detuning eo ^ — l/2/ii^*a. 
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Fig. 1. Diagrammatic series contributing to the atom-dimer 
T-matrix and a schematic representation of the Skorniakov- 
Ter-Martirosian integral equation (12). External propagators 
are included to guide the eye, they do not form part of the 
T-matrix. Straight and wavy lines denote atomic and dimer 
propagators, respectively. 



3 Atom-dimer scattering 

Knowledge of atom-dimer interaction parameters is neces- 
sary for the correct description of an atom-molecule mix- 
ture on the BEC side of the Feshbach resonance. The 
momentum-space formalism for the three-body problem 
with short-range interactions was first demonstrated in 
the calculation of the neutron-deuteron scattering length 
(total spin S = 3/2 in this case corresponds to our t- 

scattering problem) [26]. The coordinate formulation 
can be found in Ref. [18] where the atom-dimer scattering 
length was obtained in the mass-imbalanced case. Here we 
extend these results to higher partials waves, finite colli- 
sion energies, and finite Feshbach resonance width. 

Let us denote the atom-dimer scattering T-matrix by 
T(k, /co; p,po), the arguments of which imply that the in- 
coming four-momenta of the atom and the molecule are 
(k, ko) and (— k, — /cq), and the outgoing ones are (p,Po) 
and (— p, E — po), respectively. In Fig. 1 we show the di- 
agrammatic series for T, the summation of which results 
in the Skorniakov-Ter-Martirosian integral equation [26] 
(see also Ref. [27]) 

T(k,A:o;p,Po) = -g^Z G^^-k - E - ko - po) 

^^-^G^(q, go)G;(-p -q,E-po- qo) 

xD(-q,T;-go)T(k,A:o;q,go). (12) 

Equation (12) is formally identical to the equal-mass wide- 
resonance one, the difference being hidden in the propaga- 
tors and the factor Z, which serves for correct normaliza- 
tion of external propagators (see Appendix A). The atom- 
dimer elastic scattering amplitude is proportional to the 
on-shell T-matrix: 



/(k,k') = -^T(k,fcV2mt;k',fcV2mt), 



27r 



(13) 



where fis = Mm^/{M + m^) is the reduced mass of the 
atom-dimer system, and A: = |k| = |k^|. Hereafter /, /g. 
Si, and ai refer to the atom-dimer scattering parameters, 
the two-atom interaction being described by a and i?*. 

Integration over qq in Eq. (12) may be carried out by 
closing the complex contour in the lower half plane. The 
integration picks up the contribution from the simple pole 
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^K-^Li mixture Equal masses ^^K-^Li mixture 




Fig. 2. (color online). Atom-dimer s, p, and d-wave scatter- Fig. 3. (color online). Partial K-KLi atom-dimer cross-sections 
ing phase shifts vs. £^coii/|eo|. Solid, dashed, dot-dashed, and ^r, (sohd), (dashed), and (Jd (dotted) in units of vs. col- 
dotted lines correspond to W/a = 0, 1/16, 1/4, and W = a, re- Vision energy for different detunings R^a. 
spectively. In the homonuclear case we show 5d only for i?* = 0. 



of at qo = /2m^. The scattering phase shifts are on- 
shell quantities and we let = k'^/2m^, po = /2m^^ 
and the total energy E — /c^/2/i3 + eo- The remaining on- 
shell condition, |p| = |k|, should be implemented at the 
end of the calculations. 

The kernel of the resulting three-dimensional integral 
equation has a simple pole at |q| = |k| hidden in the dimer 
propagator. We make it explicit by defining functions / 
and h\ 



/(k,q) 

j2 _ J^2 _ 



D{ci,E-qy2m^l^f^ q 



2mi 



h{k, q) = (k^ - q^)D{q, E - q^/2m^)/Aii. 



(14) 



Both / and h are not singular at |q| = |k|, and / is cho- 
sen such that /(k,kO = /(k,kO, for |k| = |k'|. Finally, 
we note that Eq. (12) conserves angular momentum and, 
therefore, can be written as a set of decoupled equations 
for each partial wave 



where we define 



(15) 



9e{k,p) 



1 



fcV2m|-pV2"^t), 
(16) 

where x is the cosine of the angle between k and p. Partial 
atom-dimer scattering amplitudes are related to solutions 
of Eq. (15) by the equation fi{k) = fi{k^ /c), and the cor- 
responding phase shifts 5^ are deduced from Eq. (2). 

In Fig. 2 we plot the 5-, p-, and d-wave phase shifts as 
functions of the collision energy £^coii = ^^/2/i3 for differ- 
ent detunings R* /a. We write the phase shifts as 5q = 5s^ 



8i = 5p ^ and 82 = Sd- The results are shown for two 
mass ratios: m^/m^ = 6.64 (left) and m^/m^ = 1 (right). 
We keep the same vertical scale in both graphs, and one 
can see that the atom-dimer interaction in the heteronuc- 
lear case is stronger in every considered channel. Looking 
at the low-energy asymptotes of the phase shifts in the 
wide resonance case (i?* = 0) we see that passing from 
m^/m^ = 1 to m^/m^ = 6.64 the atom-dimer 5- wave scat- 
tering length increases from a^^ ~ 1.18a to a^^ ~ 1.98a, 
consistent with Refs. [26] and [18]. At the same time the p- 
wave scattering volume increases by more than an order of 
magnitude from v^^ ~ — 0.95a^ to v^^ ~ — lO.la^, which is 
apparently due to the vicinity of the resonance at the crit- 
ical mass ratio m^jm^ ~ 8.2 [19]. Although the mass ratio 
for the K-Li case is quite a bit smaller, our results indi- 
cate that for sufficiently small detuning one has a strongly 
marked p-wave K-KLi scattering resonance. Indeed, for 
i?* = the p-wave phase shift reaches the unitarity value 
7r/2 at a relatively small collision energy £^eoii ^ 0.1|eo|. 

In Fig. 2 we also see that the atom-dimer interaction 
decreases with detuning. We attribute this to the fact 
that at larger i?*/a the light atom spends more time in 
the closed-channel molecular state, and consequently con- 
tributes less to the atom-dimer exchange interaction. In 
a sense, increasing i?*/a is similar to increasing the mass 
of the light atom (decreasing the mass ratio): the heavier 
the atom, the weaker the exchange interaction. 

Although the p-wave resonance becomes less pronoun- 
ced near a narrow resonance, in a K-Li mixture the p-wave 
atom-dimer interaction can be strong, which is demon- 
strated in Fig. 3, where we plot the partial wave cross- 
sections (Ji(k^ = 47r(2-^ + l)/c~^ sin^ ^^(/c). We clearly see 
that for detunings i?*/a < 1 the p-wave partial cross- 
section either exceeds or is comparable to as in a wide 
range of collision energies. 

For comparison, in Fig. 4 we present partial atom- 
dimer cross-sections in the homonuclear case. We see that 
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Equal masses 




Fig. 4. (color online). Atom-dimer s-wave (black), p-wave 
(blue), and d-wave (purple) scattering cross sections for the 
homonuclear gas. Solid, dashed, dot-dashed, and dotted lines 
correspond to i?*/a = 0, 1/16, 1/4, and i?* = a, respectively. 



the 5-wave contribution is always dominant and the func- 
tional form of (Js/o? is fairly insensitive to the detuning. 

We have calculated the scattering parameters for sev- 
eral higher partial waves. Their contributions rapidly de- 
crease with I and it is worth plotting only the d-wave phase 
shifts (see Fig. 2) and scattering cross-sections (Figs. 3 and 
4). The d-wave contribution is comparable to the s- and 
p-wave ones only in the heteronuclear case and for rela- 
tively high collision energies |eo|. In Figs. 2 and 4 the 
(i-wave contribution for the homonuclear case is plotted 
only for i?* = as for finite detunings the curves are even 
closer to the horizontal axis. 

We have already discussed the atom-dimer scatter- 
ing length and scattering volume for m^jm^ ^6.64 and 
m^/m^ = 1 in the case i?* = 0. In Figs. 5 and 6 we plot 
these quantities and the effective range parameters r^d and 
/Cad [atom-dimer analogues of ro and /cq defined in Eqs. (3) 
and (4)] versus the detuning R* /a. Dotted lines in these 
graphs are obtained by using a perturbation theory in the 
limit of a very narrow resonance, ^ ^ 0, when the atom- 
dimer T-matrix can be obtained by summing the first few 
diagrams in Fig. 1. This gives an expansion in powers of 
^Ja/ ^ 1. The first two terms in the expansion of a^d 
and 'L'ad and the leading terms for the effective range pa- 
rameters r^d and k^^ read 



Ms 



-4i?*^ I 1 




--a^R^^ 
3 



a /i3 V 2/i3 



^^K-^Li mixture 

0^ 




Fig. 5. (color online). Atom-dimer scattering length ttad, s- 
wave effective range Tad , p-wave scattering volume v^^d , and the 
p-wave effective range parameter /Cad in units of corresponding 
powers of a vs. i?*/a for m^/mi = 6.64. Solid lines are exact 
and dotted lines are approximate results (17)- (20) valid in the 
limit i?* ^ a. 

Equal masses 

of 




Fig. 6. (color online). Same as in Fig. 5 but for the homonuc- 
lear case. 



3.1 Born-Oppenheimer approximation 

It is instructive to consider the enhancement of the atom- 
dimer scattering and the appearance of trimers for suf- 
ficiently large mass ratios in the Born-Oppenheimer ap- 
proximation [28]. This method was introduced in Ref. [15] 
to study Efimov physics in the system of one light and 
two heavy particles. Although the Born-Oppenheimer ap- 
proximation is not exact, it serves well to illustrate the 
essential physics leading to the resonant enhancement of 
the p-wave scattering. Here we extend it to the case of a 
resonance of finite width. 

The method takes advantage of the large mass ratio by 
assuming that the state of the light atom adiabatically ad- 
justs itself to the distance R between the heavy fermions. 
The wavefunction of the light atom can be written in the 
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-^±{R)\r-R/2\/R -^±{R)\r+R/2\/R 
M±ir)^ 1 WJWl ± 1 , o/ol • (21) 



■R/2\ 



\r + R/2\ 



It satisfies the free-particle Schrodinger equation with the 
energy 

e±{R) = -^l{R)/2m^. (22) 

The singularities of ipR^±{r) at vanishing r = r ± R/2 
satisfy the Bethe-Peierls boundary condition [29] 

[rip]~/rip I = m±{R) cotd^ [ii^±{R)] 

= -l/a + i?*4(i?), (23) 

where the light-heavy 5- wave phase shift, again denoted 
by 5 is calculated at the light-heavy collision energy 
e±{R). The equation for k,±{R) is obtained by applying 
the boundary condition (23) to the wavefunction (21): 



i^±{R) T exp [-K±{R)R] /R=l/a- 



R'^kKR). 



(24) 



The second step of the Born-Oppenheimer method con- 
sists of solving the Schrodinger equation for the heavy 
fermions by using e±{R) as the potential energy surface. 
Let us denote the corresponding heavy-fermion wavefunc- 
tion by (t){R). Since the total three-body wavefunction, 
proportional to the product (l){R)iljR^±{r)^ should be an- 
tisymmetric with respect to the permutation of the heavy 
fermions, the symmetry of depends on the choice of sign 
in Eq. (21). As '0h,+ (^) is symmetric with respect to the 
permutation R ^ —R, the heavy- atom wavefunction ^ 
is antisymmetric and describes odd atom-dimer scatter- 
ing channels. Accordingly, the lower sign in Eqs. (21-24) 
corresponds to even channels. We see how the compos- 
ite nature of the dimer leads to the ^-dependent effective 
atom-dimer potentials: by solving Eq. (24) one arrives at 
a purely attractive e+(i?) for odd channels and purely re- 
pulsive e-{R) for even ones. 

From the viewpoint of the radial Schrodinger equation 
it is convenient to introduce the total effective potential 
for each (j)i{R)'. 



V,{R) 



K {R) - e(oo) -h + 1) /mti^^ (25) 



which includes the centrifugal barrier and shifts the thresh- 
old to zero by subtracting the dimer binding energy 



e(oc) = -( vTT4R*7a - if/^m^R 



.*2 



(26) 



In the limit ^ Eq. (26) reduces to Eq. (11). 

For the p-wave atom-dimer interaction the central is- 
sue is the competition between the attractive exchange 
potential e+ oc 1/ and the repulsive centrifugal barrier, 
which is inversely proportional to m^. In Fig. 7 we show 
Vi {R) in the limit of vanishing detuning for different mass 
ratios. Remarkably, for m^/mi ^ mj^/rrii^i this potential, 
being repulsive in both limits R <^ a and R^ develops 
a well at distances R ^ a. For m^/m^ > 8.2 the depth 
of this well is enough to accomodate a trimer state with 
unit angular momentum [19] and for somewhat smaller 




Fig. 7. The Born-Oppenheimer atom-dimer effective potential 
Vi{R) [in units of l/2m4.a^] in the wide resonance case (R* = 
0) for mass ratios m^/mi = 5 (dashed), mj^/m-^-^ (solid), 8.2 
(dash-dotted), and 13.6 (dotted). 




Fig. 8. The Born-Oppenheimer atom-dimer effective potential 
Vi{R) [in units of l/2m^a^] in the K-Li case for i?* = (solid), 
R* — a/16 (dashed), i?* — a/4 (dash-dotted), and R* — a 
(double-dot dashed). The dotted line is the centrifugal barrier. 



mass ratios the presence of the well leads to the resonant 
enhancement of the p-wave interaction. 

The effect of finite R* is to decrease the strength of the 
exchange potentials e±. In Fig. 8 we show Vi{R) in the K- 
Li case for different values of the detuning R"" /a. One can 
see that the p-wave attraction becomes less pronounced 
and the well eventually disappears with increasing i?*/a. 

It is important to distinguish the p-wave trimer for 
m^/mi > 8.2 from Efimov t rimers. The former exists 
only for a > and is a result of the peculiar competi- 
tion between the exchange potential and the centrifugal 
force at distances of the order of a, which determines its 
size. In contrast, the Efimov effect occurs at larger mass 
ratios, m^/m^ > 13.6, when the effective potential at dis- 
tances ^ a is no longer repulsive. Then, in the Born- 
Oppenheimer description the heavy atoms fall to the cen- 
ter in an attractive 1/i^^-potential. This is accompanied 
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Fig. 9. (color online). The integrated column density for K- 
KLi atom-dimer scattering in arbitrary units. The collision en- 
ergies are ^coii = 0.05|eo| (left) and ^coii = 0.25|eo| (right). In 
both cases i?* = a/4. Backward direction corresponds to neg- 
ative z. For presentation purposes we imitate a small thermal 
smear. 



by the appearance of an infinite set of Efimov states, irre- 
spective of the sign of a. 

Atom-dimer scattering in even channels is described 
by the potential e+(i?), which is defined at distances R > 
a. It has a form of a purely repulsive soft-core potential, 
which increases with the mass ratio and decreases with i?* , 
consistent with the exact results above on s- and d-wave 
atom-dimer scattering. 



3.2 Interference of s- and p-waves 

Now we would like to discuss one of the implications of 
the channel dependent atom-dimer interaction, a pecu- 
liarity, which is strongly pronounced in the K-Li mixture. 
In this case the p- and s-wave phase shifts are comparable 
in magnitude and can be large (see Fig. 3). It is thus not 
necessary to go to very high collision energies for observ- 
ing the quantum interference between these partial waves 
[30]. 

Let us consider a gedankenexperiment in which a cold 
thermal cloud of KLi dimers collides with a cloud of K 
atoms at collision energies below the dimer break-up thresh- 
old. The measurable quantity is then the angular distribu- 
tion of scattered dimers (or atoms), which is proportional 
to the differential cross-section. We can write it in terms 
of the phase shifts by using Eqs. (1) and (2): 



1 r • 2 c- 



cos{Sp — Ss) sin 5s sin Sp cos 
9 sin^ Sp cos^ 6>] + . . . 



(27) 



Here k is the relative atom-dimer momentum and the an- 
gle 6 is measured with respect to the collision axis, which 
we denote by z. The dots signify the contribution of higher 
partial waves. We have checked that they can be safely ig- 
nored. 

The first term on the right hand side in Eq. (27) gives 
the well-known spherically symmetric scattering halo. The 
last term corresponds to the pure p-wave scattering. It 
contributes equally to the forward (0 < < 7r/2) and 
backward (7r/2 < ^ < tt) directions, but vanishes in the di- 
rection perpendicular to z. The second (interference) term 
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Fig. 10. The contrast vs. collision energy for detuning R* /a 
equal to (solid), 1/16 (dashed), 1/4 (dot-dashed), and 1 
(double-dot dashed). Dotted line is the homonuclear case result 
for R* = 0. 



favors either backward or forward scattering. In Fig. 9 we 
simulate an absorption image (column density) of scat- 
tered particles that initially moved in the positive z direc- 
tion. In the K-Li case backward scattering dominates at 
small collision energies while forward scattering is favored 
at higher energies, when Sp — Sg > 7r/2. 

In Fig. 10 we plot the contrast, defined as the normal- 
ized difference between the numbers of particles scattered 
forward, A/'+, and backward. A/"-, as a function of collision 
energy for different detunings R* /a. For comparison we 
also present the homonuclear wide-resonance case (dotted 
line). We see that in this case backward scattering always 
dominates, the highest contrast achieved for £^^011 ^ |eo|/3. 

The collision experiment described above requires the 
ability of manipulating atoms and molecules individually, 
which points to an advantage of heteronuclear mixtures 
- in the heteronuclear case different atomic species feel 
optical potentials in a different manner, and this obviously 
holds for the two components of the corresponding atom- 
molecule mixture. 

As far as the energy scale is concerned, in the K-Li 
case with i?* = a = 100 nm the dimer binding energy 
given by Eq. (11) equals |eo| ~ 1.8 /iK, and it decreases to 
200 nK when a = 400 nm. We also mention, for reference, 
that the relative atom-dimer velocities corresponding to 
the collision energy £^^011 = |eo| in these two cases equal 
3.7 cm/s and 1.2 cm/s respectively. 



4 Dimer-dimer scattering 

Dimer-dimer interaction parameters are crucial for the de- 
scription of the BCS-BEC crossover in the BEC-limit, i.e. 
when the gas of molecules is dilute. In the lowest order 
the chemical potential, condensate depletion, and speed 
of sound in the EEC of dimers are determined from the 
density and dimer-dimer scattering length a^^ in the same 
manner [31] as in the usual Bogoliubov theory of dilute 
Bose gases. 
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In the case of a homonuclear mixture near a wide res- 
onance the dimer-dimer scattering length equals a^d ~ 
0.6a. This number was obtained in Ref. [32] by solving 
an integral equation derived directly from the four-body 
Schrodinger equation in coordinate space. Later the result 
was confirmed by diagrammatic approaches [31,33] and by 
Monte-Carlo and variational techniques [34,35]. In the he- 
teronuclear case the molecule-molecule scattering length 
was also calculated in the case of a wide interspecies reso- 
nance [36,35,37]. Von Stecher et al. [35] also computed the 
dimer-dimer effective range. The inelastic scattering and 
the formation of Efimov trimers in dimer-dimer collisions 
for m^/mi > 13.6 is discussed in Ref. [38]. 

A qualitative summary of the results cited above is 
that the dimer-dimer interaction can be thought of as 
a soft-core repulsion. It strengthens with the mass ratio, 
but is always weaker than the 5- wave repulsion between a 
dimer and a heavy atom. This picture can be understood 
from the Born-Oppenheimer analysis, when one assumes 
that the wavefunction of the two light fermions is given by 
the antisymmetrized product of '0R,+ (ri) and ipR,-{r2) 
[see Eq. (21)], which is antisymmetric under the permuta- 
tion of the heavy fermions. Accordingly, the heavy-atom 
part of the wavefunction should be symmetric, consistent 
with the fact that only even scattering channels are al- 
lowed between identical bosons. The Born-Oppenheimer 
potential energy surface is given, independent of the an- 
gular momentum, by the sum e+(i?) + e-{R). It increases 
with decreasing the mass of the light atom and is repul- 
sive, but not as strong as the atom-dimer 5- wave potential 
e_(i?). 

As far as we know, higher partial waves in scatter- 
ing of bosonic dimers have not been studied, but it has 
been shown that the ground state of four t ~ t ~ i ~ i 
fermions in an anisotropic harmonic potential has zero 
angular momentum, independent of a and the mass ratio 
[35,39]. Besides, the qualitative Born-Oppenheimer anal- 
ysis does not provide arguments for any resonant enhance- 
ment of higher partial waves. We thus conjecture that the 
s-wave channel should dominate the dimer-dimer interac- 
tion, at least for sufficiently small collision energies. 

The aim of this section is to compute the dimer-dimer 
5- wave scattering length add for the homonuclear and hete- 
ronuclear cases (having in mind the Li-K mixture) taking 
into account the finite width of the Feshbach resonance. 
Our derivation follows Refs. [31,37] where the problem 
was studied in the regime of small detuning. 

We consider the scattering of two dimers with four- 
momenta (0, eo) into dimers with (±p, eo±Po) and project 
onto the s-wave (average over directions of p). The four- 
body T-matrix with these kinematics is denoted T{p,po). 
Similarly to the three-body case the four-body T-matrix 
consists of an infinite sum of diagrams, which may again 
be reduced to integral equations. In order to perform this 
summation, we first construct the sum of two-dimer ir- 
reducible diagrams beginning and ending in two dimer 
propagators. These are the diagrams that cannot be di- 
vided in two by cutting only one pair of dimer propaga- 
tors (external lines are excluded from the summation). 



Fig. 11. Diagrammatic representation of the integral equation 
(28), which relates the dimer-dimer T-matrix with the sum of 
all the two-dimer irreducible diagrams, T. 




Fig. 12. The integral equation satisfied by the sum of the 

two-dimer irreducible diagrams, in which one of the outgoing 
dimers is split into its constituent parts. 



The corresponding 5-wave averaged sum is denoted by 
r{q,qo;p,Po), where the four-momenta of the incoming 
[outgoing] dimers equal (ibq, eo ± go) [(ip? ipo)]- The 
equation for the T-matrix then reads 

T(p,Po) =/^'r(0,0;p,po) + ^ f q^dqdqoT{q,qo) 



xr{q, qo;p,po)D{q, eo + qo)D{q, eo - qo) 



(28) 



and is illustrated in Fig. 11. As in the three-body case, the 
prefactor of the first term on the right hand side serves 
for the correct normalization of external propagators. In 
order to avoid poles and branch cuts we solve Eq. (28) by 
rotating the contour of the go-integration to the imaginary 
axis [31]. 

The sum of the two-dimer irreducible diagrams is cal- 
culated as follows. We first sum the two-dimer irreducible 
diagrams which end in a dimer and two fermionic atoms. 
We denote the sum of such two-dimer irreducible diagrams 
by x(^5 ^o; Pi5 P2)5 where the two incoming dimers have 
four-momenta (zbq, eo ± go), and the outgoing t [|] atom 
may be put on-shell with four- momentum (pi,p^/2m^) 
[(P2 5^2/2^4.)]- By energy- momentum conservation the out- 
going dimer then has four- momentum (— pi — p2,2eo — 
Pi/2m^ — P2/2mi). In Fig. 12 we illustrate the integral 
equation satisfied by x- The equation itself reads 



x(^,^o;Pi,P2) 



Gi q-pi,eo +90 



j1 

2m^ 



xG^ |^-q-p2,eo - go 



Pi 



2m; 

Q + Pi +P2,2eo 



[(q,go) ^ -( 



q, go)] I 



2mi 2m^ 2m; J 



xD ( Q + pi,2eo - - 7^ ) x(g,go;Pi,Q) 



Q + pi+p2,2eo 



2m; 2m^ 

Pi pI 



xD Q + p2,2eo 



2m>^ 2m^ 2m; 

p^ \ 

^ * x(^,^o;Q,P2) 



2m^ 2m; 



(29) 



where the frequency integration in the closed loop of the 
iterated term is already performed. The configurational 



J. Levinsen and D. S. Petrov: Scattering in fermionic mixtures near a narrow Feshbach resonance 



9 




Fig. 13. r expressed in terms of x (see text). 



space of Eq. (29) is in fact three-dimensional. It consists 
of the moduh of the vectors pi and p2, and the angle 
between them. The pair (g, qq) enters parametrically. In 
order to express F in terms of x, it is advantageous to 
separate out the simplest diagram, in which the dimers 
exchange identical atoms. Then, the remaining diagrams 
in r are obtained by closing the fermionic loop in x (see 
Fig. 13). 

The relation between F and x is 




1 f d^pi d^p2 f 
-27 (2^(2^1^^ ^P-Pi,^o+Po 



A. 
2m-f 



xG-f ( p + P2,eo -po 



2m J 



xD pi +P2,2eo 



2mt 



A. 

2m\ 



[(p,Po) ^ -(p,Po)] 



x(^,^o;pi,P2), (30) 



where the factor ^ is needed for correct counting of dia- 
grams. The quantity /^^^^ is the first diagram on the right 
hand side of Fig. 13 and is given in Appendix B. 

The dimer-dimer scattering length is related to the T- 
matrix by 

M 

«dd = ^T(0,0). (31) 

Fig. 14 shows our results for the dimer-dimer scattering 
length in the equal mass case and for the Li-K mixture. 
In the limit of small detuning we recover the results [32, 
36,35,37] 

add = 0.60, m^/rui = 1, (32) 

ttdd = 0.89, m^/m^ = m^/m^,. (33) 

In the opposite limit the diagrammatic expansion becomes 
perturbative as is the case for the atom-dimer problem dis- 
cussed in Sec. 3. The dominant contribution to the dimer- 
dimer T-matrix is provided by r^^\ Including also the 
next order, we find [40] 



0.13, = 171^ 

0.23, m^/mi = 6.64 



R* > a. 
(34) 

The first term on the right hand side of Eq. (34) has been 
derived in the equal-mass case in Ref. [23]. 




5 Relaxation rates 

The weakly bound dimers that we are considering are in 
fact molecules in the highest rovibrational state. They can 



Equal masses 
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Fig. 14. The dimer-dimer scattering length vs. /a for equal 
masses (solid line) and for the Li-K mixture (dashed). The 
dotted lines correspond to the asymptote (34) valid in the limit 
R* > a. 



undergo relaxation into deep bound states in collisions 
with each other or with unbound atoms. The process is 
local as it requires at least three atoms to approach each 
other to a distance comparable to the size of the future 
molecular state, i.e. ~ Rq. The released binding energy is 
of the order of l/m^i?g and is much larger than all other 
energy scales in the problem including the height of the 
trapping potential. Thus, the relaxation products are lost. 

Although the relaxation is a short-range phenomenon, 
it can be treated in the zero-range approximation. For 
wide resonances in the Efimov case, i.e. for bosons or for 
fermions with m^/m^ > 13.6, the Efimov physics is well 
described by the motion of three atoms in an effective at- 
tractive 1/R'^ potential [41,42]. The three-body wavefunc- 
tion can be separated in an incoming wave and an outgo- 
ing one, and the relaxation process is taken into account 
by adding an imaginary part to the three-body parameter 
[43]. It fixes the ratio of the corresponding incoming and 
outgoing fluxes. The physical range of the potential R^ 
does not enter the resulting relaxation rate constant. 

The suppression of relaxation in the non-Efimovian 
cases {i.e. the ttl system of fermions with m^/mi < 
13.6) originates from the centrifugal barrier for identical 
fermions, which, in turn, leads to the repulsive effective 
three-body 1/R^ potential [14,44]. In order to recombine, 
the atoms have to tunnel under this barrier to distances 
^ i?e- The zero-range approach in this case is perturba- 
tive. It uses the unperturbed few-body wavefunction to 
predict the probability of finding three atoms at small 
distances and gives the functional dependence of the re- 
laxation rate constant on the scattering length for a given 
mass ratio [17,36]. If the relaxation rate constant is known 
for a certain a, one can predict its value for any other 
a > Re. 

To be more specific, let us demonstrate how one can 
estimate the atom-dimer relaxation rate in the case of a 
wide resonance, for example in 5- wave collisions. For an 
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atom and a molecule in a unit volume, the probability of 
finding them within the distance a from each other equals 
(we assume that there is no 5-wave atom-dimer reso- 
nance). At distances smaller than a, the three-body wave- 
function (in the center-of-mass reference frame) factorizes 
into lP'(Ri,R2,r) oc p^^~^^{Q)^ where the hyperradius is 
p = v^(Ri - R2)2 + m;/(2m; + m^)(2r - Ri - Rs)^, f2 
is a five-dimensional set of all the remaining coordinates 
( hyper angles ) , and ^ is a normalized hyperangular wave- 
function. The power Us for the Ht-system is given by the 
root of the transcendental equation [14] 



^ ' 2 sin(2(^) cos(7n/s/2) 







(35) 



in the interval —1 < Rez/g < 3. In Eq. (35) (j) is defined 
as (/) = arcsin [m^/{m^ + rni)]. For m^/mi = 1 we obtain 
Us ~ 1.166 and for m^/mi = 6.64 we get z/^ ^ 2.02. 

The probability of finding the three atoms at hyper- 
radii smaller than p scales with p as P{p) oc \p^^~^\'^p^. 
The last term is the volume factor of the six- dimensional 
configurational space of the three-body problem (in the 
center-of-mass reference frame). Thus, the total probabil- 
ity of finding the three atoms in the relaxation region is 
- a^P{Re)/P{a) ^ a^{Re/a)'^''^+^. The relaxation rate 
constant is obtained by multiplying this probability by the 
frequency with which the relaxation process takes place 
once the atoms are within the range of the potential. It 
can be estimated as ^ 1/m^i^g. Finally, for the rate con- 
stant we obtain 



a 



ad 



Re f Re 



2vs + l 



(36) 



In the case of a narrow resonance the three atoms can 
approach the recombination region either as free atoms or 
as a closed- channel molecule and an atom. One can show, 
however, that in this case the probability of the former is 
much smaller than the probability of the latter. Indeed, 
let us consider two atoms, t and |, in the center-of-mass 
reference frame. The state of the system is given by 

Demanding that Eq. (37) be an eigenstate of the Hamilto- 
nian (7) with energy E we get two coupled equations for 
i/^k and 00, one of which in coordinate space reads 



__^^(R) + ^0o^(R) = ^V^(R). (38) 

From Eq. (38) one can see that the singularity of '0(R) at 
the origin is related to 00 by 



'0(R ^ 0) = 0o/V47ri?*i?, 



(39) 



where we have used the second of Eqs. (8) to express g in 
terms of i?* . 



Using Eq. (39) it is straightforward to show that the 
wavefunction of the weakly bound molecular state is given 
by 

0o,b = = Jl-ll^flVm-/a (40) 



and 



^|;^,{R) = Vl - Z ^ h/2tt eyiY>{nR) / R, (41) 



where = \^2peo [see Eq. (11)] and Z is defined in Eq. (65). 
We see that the probability of finding the atoms in the 
open channel equals 1 — Z and is small in the regime of 
intermediate detuning. Therefore, as i^*/a oo the re- 
laxation rate constant tends to a constant value [45] 



(42) 



which corresponds to the relaxation in collisions of atoms 
and bare molecules. 

Equation (39) can be used in a more general situa- 
tion as it applies to a pair of atoms when they are very 
close to each other. Even in a system of more than two 
atoms and/or in an external potential we can look at a 
particular pair of atoms and observe that the probability 
of finding them in the open channel at separations smaller 
than R equals \^(R')\^47rR'^dR' = \(j)Q\^R/R\ i.e. in 
the case R <^ R"^ it is much smaller than the probability of 
finding them in the closed channel. In particular, we can 
conclude that locally, when three atoms are at the hyper- 
radius p <C i^*, they can be considered as an atom and 
a closed-channel molecule. The bare interaction between 
them is neglected in Eq. (7) as it is assumed non-resonant. 
Their induced interaction (via the exchange of the open- 
channel atoms) has a Coulomb form [22,46] and can also 
be neglected at very small distances. 

In order to estimate a^^ for narrow resonances in the 
regime of small detuning, R* <C a, we should slightly mod- 
ify the speculations that lead us to Eq. (36). At distances 
p ^ R* the three-body wavefunction behaves practically 
in the same manner as in the wide resonance case. The de- 
viation is important at distances smaller than i^*, where, 
as we have just mentioned, the three-body wavefunction 
describes a non-interacting atom and a bare molecule. The 
rate constant reads 



^ ad 



(43) 



where the first factor is the relaxation rate for an atom 
and a bare molecule confined to a volume of size i^*^ and 
the second factor is the probability to find three atoms 
in this volume. We see that the ratio rjs = ^s^/^s^hare 
interpolates between rjs ^ (i^*/a)^^^+-^ for small i^*/a to 
7^5 = 1 for large R^ /a. 

The atom-molecule relaxation rate in the p-wave chan- 
nel can be estimated in the same fashion. The difference 
from the s-wave case is the additional factor (/ca)^, where k 
is the relative atom-molecule momentum. It enters due to 
the unit angular momentum when we calculate the proba- 
bility to find the atom and the molecule at distances ~ a. 
For the same reason the relaxation rate constant in the 
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Fig. 15. The diagrams leading to the first order corrections to 
the atom-dimer T-matrix, Eq. (47). 



collision of an atom and a bare molecule is now momen- 
tum dependent, i.e. 



(44) 



We also have to take into account the angular momentum 
when we calculate the relaxation rate of an atom and a 
bare dimer confined to distances R* . It now equals ~ 
^pj3are(^)/^^^*^' whcrc wc usc the fact that Cl^pj3are(^)/^^ 
is momentum independent in the ultracold limit, kRe <C 1. 
One can show that the ratio rjp = a^^ {k) / ap^^g^^^{k) should 
behave as {R* /a)'^^p~^ for small R*/a and should tend to 
1 for large i?*/a. The power Vp is given by the root of the 
equation [18] 



cot 

+ 1 2 



sin^cos[(/)(z/p + 1)] — sin(i/p0) 
{up + 1) sin^ 



= 



)COS (/)sin(7ri^p/2) 

(45) 

in the interval —1 < Re^^p < 2. For m^/m^ — 1 Eq. (45) 
gives Vp ^ 0.773 and for m^jm^ = 6.64 we get Up ^ 0.198. 

The qualitative analysis of the relaxation rates pre- 
sented above is valid in the limits of small or large de- 
tunings. However, from the practical viewpoint the most 
interesting is the crossover region, R^/a ~ 1. To calcu- 
late the inelastic rates in the general case we add to the 
Hamiltonian (7) a weak imaginary short-range interaction 
potential between the heavy atoms and bare molecules. 



HL = -i 



Q,k,p 



• 3Z\;^k • 



V 



+ t 



.p^kfit,Q-k, (46) 



where the parameters A'f and A"^ are chosen such that in 
the extreme limit i?*/a ^ 1 the corresponding relaxation 
rate constants tend to their bare values Q^^^bare ~ '^A'f 

We can now calculate the ratios rjs and rjp for finite 
atom-dimer collision energy and for arbitrary R* /a by 
treating Eq. (46) as a perturbation to the Hamiltonian 
(7). In Fig. 15 we use the unperturbed (elastic) T-matrix 
found in Sec. 3 to construct the first order digrams con- 
tributing to the inelastic correction 5T. The explicit on- 
shell expression reads 



5Ts{k) = 
STpik) = 



iZA: 



1 + lsik) 



1 — i tan 5s{k) 
1 + 7p(fc) 



1 — it&uSpik) 



= -iA;'k'r]p{k), 



(47) 
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Fig. 16. (color online). The ratios ris (left) and r]p (right) ver- 
sus collision energy for various detunings: R* / a = 4 (dotted), 
1 (solid), 1/4 (dashed), and 1/16 (dot-dashed). 




0.1 1 

Fig. 17. (color online). The ratios r]s (top) and rip (bottom) at 
zero collision energy versus i?* /a for the equal mass case (solid) 
and the Li-K mixture (dashed). We also show with dotted lines 
the corresponding power scalings in the limit of small R* /a (see 
text) . 



where we define 



TT 

2 

7rk 



7 



V / q^dq 



V 



j q'dq 



Ksik,q) 

— k'^ 
Kpjk.q) 
— k^ 



(48) 
(49) 



These integrals should be performed by taking the prin- 
cipal value, and Ki{k^q) = fi{k^q)/[l + ikfi{k^k)] is the 
solution of Eq. (15) obtained by using the principal value 
prescription. 

Results of our calculation are summarized in Figs. 16 
and 17. In Fig. 16 we present the collision energy depen- 
dence of rjs and r^p, and Fig. 17 shows their zero energy 
values versus the detuning R* /a. According to their defini- 
tion these parameters are the suppression or enhancement 
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factors of the relaxation rate constants a^^ and ap^{k) 

compared to their bare values Q^^^bare ^p%arei^)- ^^^^ 
that in the ultracold limit the 5-wave relaxation is much 
more dangerous than the p-wave one [compare Eqs. (42) 
and (44)]. Figure 17 thus suggests that the atom-dimer 
collisions in the heteronuclear case are much less prone to 
relaxation than in the homonuclear one. Indeed, in the Li- 
K case already for i?*/a = 1 the s-wave relaxation is sup- 
pressed by three orders of magnitude. We think that this 
is promising for the system's longevity even though the 
p-wave relaxation for this /a is enhanced by an order 
of magnitude. At this point we can also make the follow- 
ing observation: in order to suppress collisional relaxation 
in the homonuclear case one should make the detuning 
as small as possible, whereas in the Li-K case there ex- 
ists an optimal detuning where the inelastic atom-dimer 
collisional losses reach their minimum. The exact value 
of such optimal detuning depends on the average momen- 
tum k and on the actual values of the bare relaxation rates 

^s%.re and <bare(^)- 

5.1 Relaxation in dimer-dimer collisions 

The dimer-dimer scattering problem is more complicated 
than the atom-dimer one in many respects. In particular, 
one has to deal not only with the relaxation channel which 
requires two heavy and one light fermion to approach each 
other to short distances but also with analogous processes 
in the heavy- light- light subsystem. In fact, in the regime of 
intermediate detuning (i^*/a ^ 1) the occupation of the 
bare molecular states is dominant, and the above men- 
tioned "three-body" channels are suppressed as they re- 
quire one of the atoms to be in the open channel. Then 
the dominant decay scenario is the relaxation in collisions 
of bare molecules. In this case the perturbation is 



Q,k,p 



blbU 



'p^Q-p 



kf^Q-k, 



(50) 



and we treat it in the same fashion as (46). Equation (50) 
correctly describes the relaxation in the limit i^* /a ^ 1 if 
we set Z\^'^ = Q^bare/^- "^^^ diagrams leading to the lowest 
order correction to the dimer-dimer T-matrix are depicted 
in Fig. 18. The on-shell ST at vanishing collision energy 
equals 



^r(o) 



2 /Vdd 



1 + 



eo + go) 



X D{-q,eo-qo)T{q,qo) 



dd^dd 
S 5 



(51) 



where T(g, go) is the solution of the dimer-dimer integral 
equation (28), and 77^^^ is the ratio of the actual relax- 
ation rate constant at a given i?*/a to its bare value (at 
R*/a = 00). This parameter is plotted in Fig. 19 versus 
the detuning R* /a and we see that this "four- atom" re- 
laxation mechanism is suppressed for small detunings. The 



power law dependence in this limit can be understood in 
a similar fashion as in the atom-dimer case. We start with 
two molecules in a unit volume. The probability to find 
them within the distance ~ a equals a^. Then, at shorter 
distances the four-body wavefunction is proportional to 
pi^4body-i^ where p is the four-body hyperradius. Thus, 
given that the four atoms are confined to the volume ~ a^, 
the probability to find them at hyperradii p < i?* equals 
(i?*/a)^^4'^°dy+'^, where we take into account that the four- 
body configurational volume in the center-of-mass frame 
scales as p^. At hyperradii p < we deal with two bare 
molecules and the relaxation rate is given by a^^^^/R*^. 
Finally, for the dimer-dimer relaxation rate constant asso- 
ciated with the "four atom" relaxation mechanism in the 
limit i^*/a ^ 1 we obtain 



(52) 



The power i^4body "lay be obtained by calculating the 
ground state energy of the four-body system in the unitar- 
ity limit in a harmonic potential which has been carried 
out in Ref. [47] for various mass ratios. We cite the follow- 
ing results [47]: z^4body ^ 0.0 for m^/m^ = 1, z^4body ~ 0.3 
for m^/rri]^ = 4, and z^4body ~ 0.5 for m^/m^ = 8. Our nu- 
merical calculations for the dimer-dimer relaxation in the 
limit i?*/a ^ 1 are consistent with Eq. (52) with these 
powers. 

Figure 19 only shows the result corresponding to the 
"four-body" mechanism of dimer-dimer relaxation (domi- 
nant at intermediate detunings) and thus presents a lower 
bound for the relaxation rate constant. We see that for 
i?*/a ~ 1 the dimer-dimer relaxation is less suppressed 
compared to the atom-dimer case. Consequently, this re- 
gion of detunings is more suitable for studies of atom- 
dimer mixtures with low concentration of molecules. Oth- 
erwise, if one wants to study mixtures with higher molecu- 
lar concentrations, it is necessary to decrease the detuning 
as much as possible. 

In the regime of small detuning the "three-body" mech- 
anisms of dimer-dimer relaxation can be as important 
as the "four-body" one. For wide resonances the "three- 
body" cases are described in detail in Ref. [17,36]. These 
results can be easily generalized to the regime of small 
detunings near a narrow resonance by counting probabil- 
ities as we did earlier in this section in order to estimate 
the atom-dimer relaxation rates. Let us present the final 
results. The dimer-dimer relaxation rate constant origi- 
nating from the 5-wave atom-dimer relaxation mechanism 
equals (up to a numerical prefactor) the one for the atom- 
dimer collisions given by Eq. (43). Namely, 



«f ad ~ afbare(^7«) 



ad 



(53) 



The dimer-dimer relaxation rate constant originating from 
the 2^-wave atom-dimer relaxation mechanism equals 



<ad~<bare(l/«)(^7«)'''''-' 

~ {Rl/m^a''){R*/af''^-\ 



(54) 



Equations (53) and (54) can be modified to describe the 
light-light-heavy relaxation mechanism. In this case the 
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Fig. 18. The diagrams leading to the first order correction to 
the dimer-dimer T-matrix, Eq. (51). 




Fig. 19. The relaxation rate in dimer-dimer collisions. 

bare rates o^^^are ^p*bare(^) should be taken for the 
collisions of light atoms and bare molecules, and the pow- 
ers Us and Up should correspond to the light-light-heavy 
three-body system with zero and unit angular momenta 
respectively. For m^/mi = rriK/mu they equal Ug ~ 1-01 
and Up ^ 0.945. 



6 Concluding remarks 

Concluding the paper, we would like to discuss several is- 
sues, which, from our viewpoint, are important for future 
studies. We believe that the p-wave atom-dimer resonance 
should significantly influence the behavior of a mixture of 
weakly bound K-Li molecules and K atoms. The "sim- 
plest" many-body problem in which the effect of the res- 
onance should be visible is the problem of a single atom 
immersed in a BEC of molecules or a single molecule im- 
mersed in a Fermi sea of atoms. Let us consider, for ex- 
ample, the former case. In the limit 'L'ad^ ^ 1? where n 
is the BEC density, we compute the correction to the dis- 
persion curve of the polaron due to the p-wave interaction 
with molecules. By using the standard diagrammatic tech- 
niques, the dispersion is found to be 



2m- 



2/i3 



(55) 



We see that in the case of a negative and large p-wave 
scattering volume t'ad, the effective mass of the polaron 



increases, which can be measured in an experiment sim- 
ilar to the one performed recently at ENS on collective 
modes of a homonuclear ^Li spin-mixture [48]. Most in- 
teresting seems to be the regime I'^adl^ ~ 1, where the 
validity of the mean-field approach is not guaranteed even 
on the qualitative level and a non-perturbative analysis is 
required. Our few-body results suggest that the crossover 
phase diagram in the K-Li case should be less "mean- 
field friendly" than the one in the homonuclear case. The 
same holds for the case of a single heavy-light molecule 
immersed in a Fermi sea of the heavy species. Variational 
analysis [49] suggests that such molecules can form a BEC 
with finite momentum (an analog of the FFLO phase) . In 
a related work [50], it was demonstrated that the phase 
diagram of a two-species bosonic mixture near a p-wave 
interspecies Feshbach resonance is quite rich and features, 
in particular, a finite- momentum superfiuid phase. 

In this paper we have shown that the atom-dimer scat- 
tering is characterized by a very strong dependence on the 
angular momentum, the effect being remarkably pronoun- 
ced for the K-Li mass ratio. To observe the interference 
between the s- and p- waves we have proposed a scattering 
experiment in which one collides a cloud of atoms with 
a cloud of molecules. Depending on the collision energy 
the scattering goes predominantly in the forward or back- 
ward directions. The contrast in this experiment is lim- 
ited to approximately 75% due to the fact that only two 
partial waves are involved. We believe that a similar ex- 
periment in a quasi- ID geometry would show a complete 
destructive or constructive interference, i.e. transmission- 
less or refiectionless ID scattering at certain values of the 
collision energy. The control over the atom-dimer scat- 
tering properties provided by the external quasi- ID con- 
finement opens up interesting perspectives for observing 
rather exotic effects. For example, in the case of a van- 
ishing transmission coefficient the propagation of a mo- 
lecule immersed in a quasi- ID gas of heavy fermions is 
suppressed, leading to degeneracies in its excitation spec- 
trum and an unusual dynamical correlation function [51]. 
Besides, as in the quasi-2D case [20], the quasi-lD con- 
finement can push the trimer state below the atom-dimer 
threshold, i.e. make it bound. The quasi-lD confinement 
of only one of the species leads to similar effects and can 
be used to control the atom-dimer scattering amplitude 
[52,53]. 

In the present paper we focus on the ^Li-^^K mixture in 
which the Feshbach molecule is bosonic and the enhanced 
p-wave interaction occurs between a boson and a fermion. 
On the other hand, our approach applies equally well to 
the case of light bosonic atoms. If, for example, the light 
particle is chosen to be ^Li, the enhanced interaction takes 
place between a fermionic atom and a fermionic dimer. In 
this case, however, relaxation in molecule-molecule colli- 
sions is not suppressed by the Pauli exclusion principle 
as two bosons and one fermion can easily approach each 
other. Nevertheless, the fact that the molecules are iden- 
tical fermions is advantageous for the system's longevity. 
Such a mixture of fermionic atoms and molecules can be 
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Fig. 20. The propagator of dimers consists of an admixture of 
closed channel Feshbach molecules (thin wave lines) and open 
channel fermionic loops. 



seriously considered for studies of the BCS-BEC crossover 
with resonant j9-wave interspecies interactions. 
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A Two-body scattering amplitude and 
propagator of dimers 

In this appendix we obtain the dressed dimer propagator 
and relate the bare parameters of the Hamiltonian, Eq. 
(7), to the physical observables of the scattering process 
between a heavy and a light fermionic atom. The propa- 
gator of the bare (closed-channel) molecule is 



Po-p'^/2M -cjo^iO' 



(56) 



The dressed propagator is obtained by resumming the di- 
agrams of Fig. 20, resulting in 



The fermion loop 77 is 



n{p,po) 



(2^)-^ 



Gt{p^Q)G^{-Q) 



(57) 



(58) 



= + ^V-Po^p'm - iO. (59) 

The cut-off at large momenta, yl, is the inverse van der 
Waals range of the potential. We have chosen to include 
an extra factor g'^ in the dressed dimer propagator. This 
is purely for bookkeeping purposes, as this propagator al- 
ways appears along with a factor of g'^ in Feynman di- 
agrams, providing an easy manner of keeping tracks of 
powers of the coupling constant. The dressed propagator 
is 



D{p,Po) = g'^ 



Po 



2M 



g'^^iA 




Fig. 21. The relationship between the scattering amplitude 
for elastic heavy-light atom scattering and the dressed dimer 
propagator. 



p^ 

2M 



iO 



(60) 



Let us now relate the bare parameters, g and uq, to the 
physical parameters a and i?* of Eq. (3). To this end we 
calculate the amplitude of elastic heavy-light atom scat- 
tering in the model Eq. (7). Let the heavy atom have four- 
momentum (k, k'^/2m^) and the light have (— k, P/2m^). 
The scattering is depicted in Fig. 21. The summation of 
diagrams takes the same form as above, Eq. (58), and 
recalling that g'^ has been absorbed into the dimer prop- 
agator, the scattering amplitude is 



T2(k,-k) 



D{0,k'^/2fi) 
27r 



/i 271 



~|~ o o ~\~ zk 

fjL g 



(61) 



Comparing with the scattering amplitude at low momenta, 
I: So-', 

(62) 

the scattering length and effective range are found to be 
.2 1 



a = 



27r 



■ ^0 



(63) 



In terms of the physical observables, the dimer propagator 
(60) then takes the form (10). 

The binding energy of the heteronuclear molecule is 
found as the energy pole of the propagator and is 



yrT4R*7^- 1 



8/ii? 



'*2 



(64) 



Since i?* > there is only one such pole on the physical 
sheet. The residue at the energy pole is needed for proper 
renormalization of the three- and four-body T-matrices, 
and is 



g'z 



(65) 



Each external dimer propagator acquires a factor ^/Z in 
the few-body scattering problems. 
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B The kernel of the 4-body integral equation 

is the result of calculating the Born diagram, 

x<.,(Q-M.| + Q.-f-|) 
p _ q £0 

2 2' 2 

j3n 



<3o + 



xG^ (-Q 



) 



2 
(66) 



with 



A 
B 
C 



2/i 8/i 8/i 

Q P , Q P 



pq , pq 

4mt 



4m^ 



= ^0 - 



2m^ 
Q q 



2m 



2m; 

Q q 

2mi 



Q q 

2/i ' 

Q P 

2/i 



(67) 
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